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Abstract. We present an event- by-event hydrodynamical framework which takes 
into account the initial density fluctuations arising from a Monte Carlo Glauber model. 
The elliptic flow is calculated with the event plane method and a one-to-one comparison 
with the measured event plane «2 is made. Both the centrality- and pr-dependence of 
the vi are remarkably well reproduced. We also find that the participant plane is a 
quite good approximation for the event plane. 



1. Introduction 

Hydrodynamical models using smooth initial conditions cannot reproduce the centrality 
dependence of the elliptic flow coefficients t> 2 , measured using the event plane method 
(see e.g. Fig. 7.5 in [T]). To study this more carefully we present here an event-by-event 
hydrodynamical framework introduced in Ref. [2] and make a one-to-one comparison 
with the experimental event plane results for the elliptic flow. Also the t>2 with respect 
to the participant plane is considered. 

2. Event-by-event hydrodynamics framework 

The initial state is obtained from a Monte Carlo Glauber model. First, nucleons are 
randomly distributed into the nuclei using a standard Woods-Saxon potential. Then 
the impact parameter b is sampled from a distribution dN/db ~ b. A nucleon i and a 
nucleon j from different nuclei collide if their transverse locations are close enough, 



where a^^ is the inelastic nucleon- nucleon cross section. For y / sJvjv = 200 GeV we 
use a^N = 42 mb. The impact parameter defines the reaction plane (RP) and the 
participant plane (PP) is obtained from the participant configuration by maximizing 
the eccentricity. Centrality classes are defined using the number of participants iV part 
and slicing their distribution into Ap art intervals as shown in Ref. [2]. 
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Figure 1. The transverse momentum spectra of positively charged pions in Au+Au 




The initial energy density e in the transverse plane is obtained by distributing 
energy around the participants using a 2-dimensional Gaussian smearing, 



where K is a overall normalization factor and a is a free parameter, which controls the 
width of the Gaussian. The overall normalization as well as the initial time r = 0.17 fm 
are fixed from the EKRT model [3]. 

For each event, we solve the ideal hydrodynamical equations d^T^ = assuming 
longitudinal boost-invariance and zero net-baryon density. We further need to specify 
an equation of state (EoS) P = P(e) to close the set of equations. Our choice is the 
EoS from Ref. [I]. The freeze-out is assumed to happen at a constant temperature 
Tp = 160 MeV, and the thermal spectra are calculated using the Cooper-Frye method. 
Hadrons are sampled from the calculated spectra and they are given to PYTHIA 6.4 
[5], which does all the strong and electromagnetic resonance decays. 

In order to make a one-to-one comparison with the experiments, the elliptic flow 
is calculated with the event plane (EP) method. The event flow vector for the second 
harmonic is 



where we sum over all particles, and the event plane angle ip2 is obtained from this as 





Q2 = XK^cos(2</) i ),p Ti sin(20 i )) 



(3) 



i>2 



&rct&n(Q 2 ,y/Q2,x) 
2 



(4) 
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Since we have only a finite number number of particles in each event, the event plane 
fluctuates around the true event plane, which would be the one determined from 
infinitely many particles. The correction TZ from these fluctuations is calculated with 
the 2-subevent method [6]. The final event plane elliptic flow is thus 

t; 2 {EP}= (cos(2(&--0 2 )))/ft. (5) 
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Figure 2. The elliptic flow of charged hadrons in Au+Au collisions at y / ijv r /v = 
200 GeV. The event- by- event calculations with respect to different reference planes 
and the calculations with smooth initial states are from [2]. The data are from the 
PHENIX and STAR Collaborations EJ LTO]. 



3. Results 

For each centrality class we have made 500 hydro runs and to increase the event statistics 
we have sampled each thermal hadron spectra 20 times. The smooth initial states are 
obtained by taking an average over 20 000 initial energy density profiles. The obtained 
transverse momentum spectra of pions are plotted in Fig. Q] The density fluctuations 
increase the number of high-pr particles since the initial pressure gradients are larger. 
However, if the fluctuation size parameter a is sufficiently large, the difference to the 
smooth result is small. 

In Fig. [2] we have plotted the elliptic flow with respect to the different reference 
planes. The EP results agree very well with the measured data at all centralities shown 
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Figure 3. The correlation of the event plane with the participant plane and the 
reaction plane in Au+Au collisions at ^/snn = 200 GcV. From [5] . 

in the figure. If we use PP instead of EP, the elliptic flow is a little bit smaller, but PP 
seems to be a quite good approximation for the EP. The ^{RP} is generally smaller 
than f 2 {EP}, but it is very close to the smooth result. Thus the fluctuations alone do 
not generate more elliptic flow, but the reference plane definition is very important. 

In Fig. [3] is shown the correlation of the EP with the PP, and with the RP. Note 
that the trivial fluctuations of ip2 around the true EP are included in the figure. The 
correlation between the EP and the PP is, as expected, stronger than that between the 
EP and the RP. 
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